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ABSTRACT 

Q^l We present two-dimensional stellar and gaseous kinematics of the inner 120 x 250 pc^ 

I . of the Liner/Seyfert 1 galaxy M 81, from optical spectra obtained with the GMOS inte- 

O ' gral field spectrograph on the Gemini North telescope at a spatial resolution of « 10 pc. 

The stellar velocity field shows circular rotation and overall is very similar to the pub- 
^ ■ lished large scale velocity field, but deviations are observed close to the minor axis 

' which can be attributed to stellar motions possibly associated to a nuclear bar. The 

stellar velocity dispersion of the bulge is 162± 15kms^^, in good agreement with pre- 
vious measurements and leading to a black hole mass of Msff = 5.5l^;^ xIO^Mq based 
I on the Mbh-o' relationship. The gas kinematics is dominated by non-circular motions 

. and the subtraction of the stellar velocity field reveals blueshifts of « — lOOkms"^ on 

the far side of the galaxy and a few redshifts on the near side. These characteristics 
can be interpreted in terms of streaming towards the center if the gas is in the plane. 
On the basis of the observed velocities and geometry of the flow, we estimate a mass 
' inflow rate in ionized gas of Ri4.0 x 10~^ Mq year~^, which is of the order of the ac- 

cretion rate necessary to power the LINER nucleus of M81. We have also applied 
the technique of Principal Component Analysis (PGA) to our data, which reveals the 
presence of a rotating nuclear gas disk within w 50 pc from the nucleus and a com- 
pact outflow, approximately perpendicular to the disk. The PGA combined with the 
observed gas velocity field shows that the nuclear disk is being fed by gas circulating 
. , in the galaxy plane. The presence of the outflow is supported by a compact jet seen 

. in radio observations at a similar orientation, as well as by an enhancement of the 

[01] /Ha line ratio, probably resulting from shock excitation of the circumnuclear gas 
by the radio jet. With these observations we are thus resolving both the feeding - via 
the nuclear disk and observed gas inflow, and the feedback - via the outflow, around 
the low- luminosity active nucleus of MSI. 

Key words: Galaxies: individual (M81) - Galaxies: active - Galaxies: Seyfert - 
Galaxies: nuclei - Galaxies: kinematics - Galaxies: jets 



1 INTRODUCTION 

It is widely accepted that the radiation emitted by an active 
galactic nucleus (AGN) is a result of accretion onto the cen- 
tral supermassive black hole (hereafter SMBH). However, 
the exact nature of the mechanisms responsible for trans- 
fering mass from galactic scales down to nuclear scales to 
feed the SMBH is st ill an open question. Theoretical stud- 
ies and simulations (IShlosman et aL II1990I : lEmsellem et al.l 
I2OO3I : iKnapenllioOsI : lEmsellem et al.ll2006l ) have shown that 



non-axisymmetric potentials efficiently promote gas inflow 
towards the inner regions jEnglmaier fc Shlosmanl liooi). 
Imaging studies have revealed that structures such as small- 
scale disks or nuclear bars and associated spiral arms are 
frequently obse rved in the inner kiloparse c of active galaxies 
(lErwin fc Spar kc 1999; Poggc & Martinil [200^ : iLaine et al.1 
I2OO3!), although Shaw ct al. (1993) raises important points 
about the identification of nuclear bars. 

While bars can be effective in transporting gas into 
the inner few hundred parsecs, the fundamental problem 



2 A. Schnorr Miiller et al. 



of how gas gets from the re down to the SMBH h as re- 
mained unsolved. Recentlv lSimoes Lopes et al.l (|2007l ) have 
shown a strong correlation between the presence of nuclear 
dust stru ctures (filaments, spi r als an d disks) and activity in 
galaxies. iGarcia- BuriUo et al.] (|2005h have argued that stel- 
lar gravity torques are also a mechanism that can drive gas 
inwards at least to a few hundred parsecs of the nuclei. 

Inward gas streaming on galactic scales has only 
been observed in a few objects. lAdler fc~Westpfah]| () 19961 ) 
found gas inflow s along the spirals arms of M 81 and 
iMundell fc Shonel (|l999l) found inflows along the bar of 
NGC4151. An important recent development has been 
the recognition of streaming motions on nuclear scales 
in both the ionized and molecular gas. In the ion- 
ized gas, inflows ha ve been ob s erved in the central 
region of NGC 1097 (iFathi et al.1 120061 ') and NGG6951 
lIStorchi-Bergmann et al 200?!') and in molec ular hydro- 
:en (Ha ) in NGC 1068 dSanchez et al] |2009|) , NGC 4051 
Riffel et all |2008|). NGC 1097 llDayies et all |2009| ) and 



Mrkl066 (iRiffel fc Storchi-Bergmannll201lD . 

With the goal of looking for more cases of inward 
streaming motions we began a project to map the gaseous 
kinematics around nearby AGN. In the present work, 
we present results obtained from integral field spectro- 
scopic observations of the nuclear region of M81, a spi- 
ral galax y with Hubb l e typ e SA(s)ab. At a distance of 
3.5 Mpc (jPaturel et al.l I2OO2I ) . corresponding to a scale of 
17pcarcsec~^, M81 harbors the nearest directly observ- 
able low-lumino si ty AGN, classified as LINER / Sey fert 1 
l|Heckmanl Il980l : iPeimbert fc Torres-PeimbertI Il98ll ). Be- 
cause of M81 proximity, the AGN has been the sub- 
ject of ma ny studies. It was found to vary w ith time 
at optical (iBower et all Il996l). radio (IHoI 119991 ) and X- 
ray wavelengths l lvomoto fc Makishimal 1200 ll ). The nu- 
clear spectrum presents broad op t ical a nd UV emission 
lines (|Peimbert fc Torres-PeimbertI 198ll) and a f e ature - 
less UV continuum ||Ho et al.lll996h . iDevereux et al.l (119971 ) 



detected a compact o ptical / UV continuum source, and 
iBietenholz et al.l (I2OO0I ) detected an ultra compact radio 
source ass ociated with a r adio j et. Its accretion physics was 
studied bv lMarkoff et al.1 l|200a ). 

The present paper is organized as follows. In Section [2] 
we describe the observations and reductions. In SectionOwe 
present the results. In section |4] we present the stellar and 
gaseous kinematics. In section [5] we present the gas flux dis- 
tributions and ratios, in section |B] we discuss our results and 
an alternative method for the analysis of the gas kinemat- 
ics using the technique of Principal Component Analysis. In 
Section 16.31 we present an estimate of the mass inflow rate 
and in Section [7] we present our conclusions. 



2 OBSERVATIONS AND REDUCTIONS 

The observations were obtained with the Integral Field Unit 
of the Gemini Multi Object Spectrograh (GMOS-IFU) at 
the Gemini North telescope on the night of December 31, 
2006 (Gemini project GN-2006B-Q-94). The observations 
consisted of three adjacent IFU flelds (covering 7 x 5arcsec^ 
each) resulting in a total angular coverage of 7 x 15 arcsec^ 
around the nucleus, with the longest extent of the fiel d cho- 
sen to be along the galaxy major axis as reported by iGoadI 



(|l976l ) (position angle PA =150°). Three exposures of 191 
seconds were obtained for each field, slightly shifted in or- 
der to correct for detector defects after combination of the 
frames. The seeing during the observation was 0.6 ", as mea- 
sured from the FWHM of a spatial profile of the galaxy ex- 
tracted along the broad Ha profile under the assumption 
that the broad-line region is unresolved. This corresponds 
to a spatial resolution at the galaxy of 10.2 pc. 

The selected wavelength range was 5600-7000 A, 
in order to cover the Ha-f [N II] AA6548,6583 and 
[S II] AA6716,6731 emission lines, observed with the 
grating CMOS R400-G5305 (set to central wavelength 
A 6300 A) at a spectral resolution of R ~ 2000. 

The data reduction was performed using specific tasks 
developed for CMOS data in the gemmi.gmos package as 
well as generic tasks in The reduction process com- 

prised bias subtraction, flat-fielding, trimming, wavelength 
calibration, sky subtraction, relative flux calibration, build- 
ing of the data cubes at a sampling of 0.1" x 0.1", and fl- 
nally the alignment and combination of the 9 data cubes. 
Unfortunately, in many of the spectra the position of the 
[S II] AA6716,6731 emission lines on the detector coincided 
with the position of one of the gaps between the CCDs, 
thus these lines could not be measured at several locations, 
as is the case of the spectrum extracted from the position 
marked as A on Fig.[T] As our data cubes are somewhat 
oversampled (relative to the seeing), a Gaussian filter with 
FWHM = 3 pixels was applied to the centroid velocity, ve- 
locity dispersion, flux and line ratio maps in order to reduce 
the noise in the flgures. 



3 RESULTS 

In Fig.[T]we show in t he upper left panel a large scale image 
of M 81 in the R band (|Cheng et al.lll997l ). In the upper right 
panel we show an image of the inner 22" x 22" of the galaxy 
obtained with the WFPC2 (Wide Field Planetary Camera 
2) through the filter F547M aboard the HST. In the mid- 
dle left panel we prese nt a structure map of the WFPC2 
HST image of M81 ( see ISimdes Lopes et al] (|2007l )). where 
a large nuclear spiral structure, as delineated by the con- 
trast provided by the dust lanes, can be seen. The rectangle 
shows the field-of-view (hereafter FOV) covered by the IFU 
observations. In the middle right panel we present an im- 
age from our IFU observations obtained by integrating the 
flux within a spectral window containing only continuum. In 
the lower panel we present three spectra of the galaxy cor- 
responding to locations marked as A, B and N in the IFU 
image and extracted within apertures of 0.1" x 0.1". 

The spectrum corresponding to the nucleus (marked 
as N in Fig.[l]) is similar to that obtained previo usly with 
the HST (|Bower et al.ll 19961 : lOevereux et al.ll2003l ). showing 
[O I] AA 6300,6363, [Nil] AA6548,6583 and [SII] AA6717,6731 
emission lines, which are broader than in the extranuclear 
spectra. The broad double peaked Ha profile reported by 
iBower et al.l (|l996l ) is not easily seen in our observations. 



^ IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Figure 1. Top left: large scale image of M81. Top right: WFPC2 image. Middle left: structure map. The rectangle shows the field of the 
IFU observation. Middle right: continuum image from the IFU spectra (flux in ergcm^s"'^ per pixel). The dashed white line indicates 
the position of the line of nodes. Bottom: spectra corresponding to the regions marked as N, A and B in the IFU image. 



suggesting that its flux has decreased since 1995, when that 
spectrum was obtained. Spectra from locations A and B 
show fainter emission in [Nil], Hq, and [SII] on top of a 
strong stellar continuum. 



4 KINEMATICS 



4.1 Stellar Kinematics 



The spectra of M81 have a high signal-to-noise ratio in the 
continuum (see Fig.[T} which has allowed us to use the stel- 
lar absortion features between 5600 A and 7000 A to derive 
the stella r kinematics. The Penalized P ixel Fitting technique 
(pPXF) (|Cappellari fc Emsellem|[20o3 ) was employed to ob- 
tain the stellar centroid velocity, velocity dispersion and the 
Gauss-Hermite moments /13 and /i4 from each spectrum. 
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Figure 2. Top left: stellar centroid velocity field (kms~^). Top 
right: stellar velocity dispersion (kms~^). Bottom left: Gauss- 
Hermitte moment h^. Bottom right: Gauss-Hermitte moment h^. 




-2 2 




-2 2 
(arcsec) 



2 2 




Figure 3. At the top, from left to right: stellar centroid velocity 
field (kms~^), model velocity field (kms~^), residuals between 
stellar and model velocity field (kms~^). Note the change of scale 
between the velocity fields and the residuals The straight white 
line indicates the position of the line of nodes. At the bottom we 
show the stellar and model velocity curve along the major axis 
(top) and the residuals along the major axis (bottom) 



using as templates the iBruzual k, CharlotI (|2003l 'l models. 
Although these models have an instrumental width (a) of 
61kms"^ at 6300A, very similar to the value of 66kms 
of our observations, we have corrected our measurements by 
this small effect, which results in a decrease of ~ 1.5 km s^^ 
in the stellar velocity dispersions. Monte Carlo simulations 
based on the best-fitting absorption spectra were carried out 
to estimate the errors in the kinematic parameters. 

The stellar centroid velocity (V,) field (top left panel 
of Fig. [21) displays a rotation pattern in which the SE side 
of the galaxy is aproaching and the NW side is receeding. 
Under the assumption that the spiral arms are trailing, it 
can be concluded that the near side of the galaxy is the 
SW, and the far side is the NE. The calculated errors in the 
velocity measurement are « lOkms"'^ over the whole FOV. 
The stellar velocity dispersion (top right panel of Fig.[21) 
decreases from decreases from ^190kms~^ at the nucleus 
to «140 kms~^ at a radial distance of ~ 5"(~ 85 pc from the 
nucleus), with an average value of 153kms~^. Calculated 
errors are «15kms~^. 

The Gauss-Hermite moments /13 and hi (bottom pan- 
els of Fig. [21 are measures of the skewness and the kur- 
tosis, respectively, of the line-of-sight velocity distribution 
(LOSVD). Errors are of the order of 0.03 for both /13 and 
/14. There is a ring- like region of «0 values at Ri20pc from 
the nucleus in the h-j. map indicating that the LOSVD in 
this region has a Gaussian profile, h-j. rises to ~ 0.05 within 



a region with radius «17pc centered on the nucleus and 
also beyond the ring, indicating a skewness of the LOSVD 
towards the red in these regions. Low values of /14 occur 
mainly in the nucleus (/i4«0) and the highest along the 
minor axis (/i4~0.06), indicating a LOSVD slightly more 
"pointy" than a Gaussian at the later locations. 

We have modelled the centro id velocity field using circu- 
lar orbits in a Plummer potential (Plummcr lQll' ) as we have 
successfully done in previous works ( Barbosa et al.l l2006l : 
iRiffel et al.ll2008l ). This model is just a kinematic model and 
does not take into account the full LOSVD of the stellar 
components. The stars are assumed to orbit in circular or- 
bits in a thin disk coplanar with the main galactic disk, in 
spite of the fact that the stellar velocity dispersion is not 
negligible. In addition we did not take into account the see- 
ing and instrumental effects as our goal is only to derive 
the mean velocity field to be subtracted from the observed 
velocity field of the gaseous component. 

The fit of the model gives a heliocentric systemic ve- 
locity of —125 ±11 kms^^ taking into account measurement 
errors in the fit, a kinematic center at horizontal and verti- 
cal coordinates of X = 0.2" and Y = QA" , respectively, rela- 
tive to the peak fiux in the continuum (adopted as the center 
of the galaxy), and a position angle for the line of nodes of 
PA = 153° ±1°. The displacement between the kinematic and 
photometric center is smaller than the seeing radius and it 
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Figure 8. Upper panel: Sample spectrum of the nucleus (black 
line) with the multi-Gaussian fit (represented by asterisks) to 
the emission lines superimposed. The individual Gaussian com- 
ponents are plotted below the spectrum along with the residu- 
als between the spectrum of the nucleus and the fitted spectrum 
(crosses). 



can thus be concluded that the two centers agree within the 
uncertainties. We have adopted a disk inclination of i = 58° , 
obtained from the aparent axial ratio (from NEeQ) under 
the assumption of a thin disk geometry, because the incli- 
nation is not well constrained in our multi-parameter fitting 
procedure. The stellar centroid velocity field, the fitted circu- 
lar model and the residuals between the stellar and modelled 
field are shown in the top panels of FigO The residuals are, 
at most locations, smaller than lOkms^^, for a total velocity 
amplitude of ~100kms~^, revealing a good fit. Neverthe- 
less, it can be noticed that the isovelocity curves of the ob- 
served velocity field show a small "s-shape" distortion when 
compared to the model, associated to the bluest and reddest 
residuals. In the bottom of Fig. [3] we show the observed and 
model velocities extracted within a pseudo-slit ~ 0.2" wide, 
oriented along the major axis at PA =153°, in order to il- 
lustrate better the fit. The systemic velocity obtained from 
the fit Vs = —125 ±11 km s~^ was adopted as the galaxy sys- 
temic velocity and subtracted from the stellar and gaseous 
centroid velocities. 



4.2 Gaseous Kinematics 

In order to measure the gas kinematics we have subtracted 
the contribution from the stellar p opulation. This was 
done using the STARLIGHT software (|Cid Fernandes et al.l 
I2OO5I ). in which the best-fitting absorption spectra is ob- 
tained through the technique of spectral synthesis. An illus- 
tration of the procedure is shown in Fig. [4] for the nucleus 
and for locations A and B (Fig.[T]). 

The gas kinematics - centroid velocities and velocity 
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dispersions - as well as the emission-line fluxes were obtained 
by the fit of Gaussians to the [Nil], Hq, [OI] and [SII] 
emission lines in the subtracted spectra. We note that the 
[S II] line could only be measured in the half of the field 
where the line profile was not affected by the gap between 
the CCDs. 

Errors were estimated from Monte Carlo simulations 
in which gaussian noise is added to the spectra. One hun- 
dred iterations were performed. Errors in the velocity, flux 
distribution and velocity dispersion measurements for the 
[Nil] and Ha emission lines have their lowest values within 
a radius of 2" and grow when approaching the borders of 
the field, where the signal-to-noise ratio is lower. The error 
values in both the velocity and velocity dispersion are re- 
spectively ~10kms~^ for the inner region and «20kms~^ 
for the outer region. For the emission line fluxes they are 
respectively « 10% for the inner region and ~ 25% for the 
outer region. In the case of the [S II] emission lines, errors 
in the velocity and velocity dispersion are «15kms~^ and 
in the fluxes « 20% over the whole FOV. For the [O I] emis- 
sion line errors in the velocity and velocity dispersion are 
~ lOkms"^ and ~ 15kms~^ respectively, while in the fluxes 
they are ~ 15%. Regions with errors in velocity higher than 
~20kms^^ were masked out. 

Within the inner 1.5" it was also necessary to fit a 
broad component to Ha. The need for this became obvi- 
ous when we examined line ratio maps such as [NII]/Ha: 
within a circular region with radius of « 1 arcsec around the 
nucleus, the ratio decreased abruptely from values of « 1.5 
outside this region to values of « 1 inside. By including a 
broad component within the inner region, the ratio between 
the narrow components [Nil] /Ha showed a smooth transi- 
tion between the two regions. In addition, by including the 
broad component in the fit, the width of the narrow compo- 
nent of Ha resulted practically the same as those of the [N II] 
emission lines, giving support to the presence of the broad 
component. As an additional "sanity check" to the above 
procedure we also verified the resulting [N II] AA6583/6548 
ratio which kept a value ~ 3.0 ± 0.2 over the whole region 
where the fit to the broad line was performed. 

An ilustration of the fit of the [N II] AA6548,6583 A lines 
and the broad and narrow Ha components is presented in 
Fig.[S] together with the residuals from the fit. Most of the 
fiux of the broad component is well reproduced by a gaussian 
with a = 862 ±37kms"^ and V = -25±60kms"^ 

Centroid velocity maps for the emission lines are shown 
in Fig.[5] and velocity dispersion maps are shown in Fig.O 
The velocity fields in the [N II] and Ha are practically identi- 
cal, and show that the gas is not in rotation along the same 
axis of the stars. If the gas were in rotation, the orienta- 
tion would be perpendicular to that of the stellar rotation. 
Since the centroid velocities for the different emission lines 
are almost identical, we chose the [N II] centroid velocity to 
represent the gaseous velocity field due to its higher signal- 
to-noise ratio over most of the field. 

In order to isolate possible streaming motions we con- 
structed a residual map by subtracting the stellar centroid 
velocity field from the gaseous one, under the assumption 
that the former represents the component of motion due 
to the galaxy gravitational potential. The residual map is 
shown in Fig.[7] along with the structure map and the [N II] 
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Figure 4. Observed spectrum (black), fitted absorption spectrum (green) and residual spectrum (blue) for three different regions. 
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Figure 5. Gaseous centroid velocity (kms^"*^) for the [Nil], Ha, [OI] and [SII] emission lines. 
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Figure 6. Velocity dispersion (kms ^) for the [Nil], Ha, [OI] and [SII] emission lines. 
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Figure 7. From left to right: structure map, gaseous centroid velocity {kms~^) and the residual between gaseous centroid velocity and 
stellar velocity field (kms~^). The vertical white line in the central panel indicates the position of the line of nodes, while the arrow in 
the right panel shows the orientation of the compact radio jet. 



velocity field previous to the subtraction of the stellar ve- 
locity field. 

The Ha and [Nil] velocity dispersion maps (hereafter 
a Ha and (tjjv//]) are very similar. In most of the field both 
o'lNii] and ana present the same values (within errors), 
usually ranging from 60kms~^ to 80kms~^. An exception 
occurs inside a radius of ~ 0.8" from the nucleus, where 
'^[Nii] is about 20kms~^ larger than ana (uncertainties 
being ~5kms~^ in this region). Also inside this radius 
""[AT//] and ana reach their highest values, «140kms~^ 
and «130kms~^ respectively. Away from the nucleus, an 
increase in the velocity dispersion is seen at two locations: 
ss5" NW and ~4" E from the nucleus, extending approxi- 
mately diagonally in the IFU towards the S. These locations 
correspond to the largest blueshifts in the gas velocity fields. 



5 LINE FLUXES AND EXCITATION OF THE 
EMITTING GAS 

In Fig.[9]we present the integrated flux distributions for the 
[Nil] A6583A, the narrow component of Ha, [OI] A6300A 
and [S II] A6716 A emission lines. While the [O I] flux distri- 
bution is compact, restricted to the inner ~ 2 arcsec around 
the nucleus, the remaining flux distributions show an elon- 
gation from NW to SE, with the Ha flux distribution sug- 
gesting the presence of spiral arms at the top and bottom 
borders of the field. In the top left of the Ha panel, to the 
north of the nucleus, there is a small region with enhanced 
emission, suggesting the presence of an HII region at this 
location. 

In Fig. [TO] we show the line ratio maps [NII]/Ha and 
[OI]/Ha obtained using the narrow component of Ha, as 
well a gas density map, obtained from the [SII] AA6717/6731 
line ratio as suming an electronic temperature of lOOOOK 
l|Osterbrock| [r989'). The [Nil] /Ha ratio is "LINER-like" - 
with [N II] stronger than Ha - over most of the IFU field. 



The only region with smaller ratios, typical of HII regions 
is observed at the location of enhanced Ha emission to 
the north of the nucleus, supporting the interpretation that 
there is an HII region there. The highest values are observed 
at the nucleus (where the [N II] /Ha ratio is « 3) and decrease 
outwards. 

The [01] /Ha ratio (central panel of Fig.llOp presents 
values typical of AGN, but there is a gradient in the values 
which range from r; 1.0 at ~ 1" SW of the nucleus to ~ 1.9 at 
^ 0.5" to the NE of the nucleus along PA « 50-60 (along the 
horizontal direction to the left in the figure). The asymmetry 
in the [O I] /Ha ratio in the nuclear region seems to be mostly 
due to an asymmetry in the [O I] flux distribution. 

The [S II] density reaches its highest value within the in- 
ner arcsecond (17pc), where it exceeds 10^ cm~^, decreasing 
to «500cm~^ at 2 arcseconds (34 pc) from the nucleus, the 
lowest measurable values are ~ 300 cm~'^ observed farther 
away. 



6 DISCUSSION 

6.1 Stellar Kinematics 

A previous study of th e stella r kinematics along PA = 157° 
by IVega Beltran et al.l (|200ll ). although on a much larger 
scale, has shown a smooth rotation pattern for the stars, 
which looks similar to ours, considering this speciflc PA. 
Nevertheless, our two-dimensional coverage allowed us to 
verify the presence of an s-shape d istortion in the iso- velocity 
contours. lElmegreen et al.l l| 19951 ) have argued for the pres- 
ence of a nuclear bar extending up to ~0.5kpc (15 arcsec) 
from the nucleus along PA « 140° , inferred from the flt of 
elliptical isophotes to a continuum J-band image of M81. 
This bar is thus oriented only 10° from the vertical direc- 
tion in Fig. [3] and extends beyond the FOV of our observa- 
tions with reaches only 7'.'5 along this PA. We may thus be 
sampling a region inside the bar and both the presence of 
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Figure 9. Maps of the [Nil], Ha, [OI] and [SII] integrated fluxes (ergcm^s^^ per pixel) 
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Figure 10. From left to right: [NIIJ/Ho? and [OI]/Ha emission line ratios and the gas density (cm"'^) derived from the [SII] line ratio. 



an s-shape distortion in the iso-velocity contours and the x- 
shaped residuals in Fig.[3]are consistent with the elongated 
pattern of stellar orbits required to support a bar. 

Considering an inclination of 58°, the amplitude of 
the rotation curve to each side of the nucleus reaches at 
most 60kms~^. This is much smaller than the velocity dis- 
persion observed over the whole IFU field, which ranges 
from 140 to 190kms~^. We thus can use the stellar ve- 
locity dispersion map (Fig. [3]) alone to obtain a virial es- 
timate for the mass of the bulge. The average value of the 
velocity dispersion inside a radius of 3.5" (the maximum 
covered by our observations along the shorter axis of the 
IFU) is 167kms~^. In order to obtain the velocity disper- 



sion of the bulge this value needs to be cor rected accord- 
ing to equation (2) of jjorgensen et aT I (1 19951) Ado pting an 
effective radius of re = 66.2" ( Baggett et al.l [l998l ) we ob- 
tain 162kms~^, which is in good agr eement with previ- 
ous measurements (.Bcifiori et al., .20091; IVega Beltran et al.l 
l200ll : lNelson fc Whittle!! 19951 ). This leads to a virial mass of 
2.1 xlO** Mq. 

The above bulge velocity dispersion can also be used to 
estimate the central b l ack hole mass using t he M-ct relation 
iTremaine et al.! 120021 : ICiiltekin et al.!!2009l ). which results 
Ms// = 5.5I2 Q xIO^Mq. This value is consistent with the 
previous virial estimate of the to tal mass, and is of the same 
order of the values obtained by iDevereux et~aLl (|2003! ) and 
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iBeifiori et al] (|2009l ') for the mass of the black hole, which are 
6.7 X lO'^ Mq and 3.7 x 10^ Mq respectively when rescaled 
to the distance we adopted in this study. 

6.2 Gaseous Kinematics and Excitation 



On kiloparsec scales, the H I gas kinematics (|Rotsl Il975l ) 
shows a similar rotation pattern to that of our stellar ve- 
locity field, and our PA = 153° of the line of nodes i s com - 
parable to the value of PA = 150° obtained bv lGoadI (|l976t ). 
Nevertheless, on the scales probed by our observations, the 
gas kinematics is very different from that of the stars. This 
was already noticed bv lGoadI (| 19761 '). who obtained the kine- 
matics for the [N II] A6583 A and Ha emission lines in the 
inner 1 kpc of the galaxy and found the gaseous velocity 
field to be complex and impossible to describe by simple 
rotation alone. A disturb ed kinematics was also found by 
IVega Beltran et all (|200ll ') in the [O III] A5007 emission-Une 
gas. But they argue that their long-slit data along PA= 157° 
suggest that there is a circumnuclear Keplerian disk of ion- 
ized gas. Our Hof kinematics along this axis does not show 
this behaviour. Probably because [O III] and H a are probing 
distinct kinematics. 

The gas isovelocity curves in our Fig.[5] suggest a kine- 
matic major axis almost perpendicular to that of the stars 
and even more than one rotation center. Nevertheless, the 
location of these apparent rotation centers, as well as the 
isovelocity contours in their vicinity, correlate with the lo- 
cations of increased velocity dispersions (Fig. [6]) suggesting 
that these regions are related to shocks in the gas. The over- 
all orientation of these features seem also to be the same as 
those of the dust lanes in the structure map. We also note 
that most blueshifts are in the far side of the galaxy, sug- 
gesting inflow towards the center if most of the gas is in 
the galaxy plane. In this case the shocks mentioned above 
could be responsible for dissipating angular momentum and 
allowing the gas to flow in towards the nucleus to feed the 
AGN. If a nuclear bar is pre s ent, a s suggested by the pho- 
tometry of lElmegreen et al.l (|l995l ') and by the distortions 
of the stellar isovelocity curves, this bar could also play a 
role in channeling gas inwards from the outer regions of the 
galaxy. This possibility is supported by the gas flux maps 
(Fig.O, which are elongated towards the bar orientation. 

If the gas is not in the galaxy plane, the observed 
blueshifts could be due to outflow to high galactic latitudes. 
Relevant to this hypot hesis is the observ a tion o f a compact 
one-sided r adio jet by Bie tenholz et al.l (|200Cll ) and, more 
recently, bv lMarkoff et al.l (12008). This jet is seen at an av- 
erage position angle of PA = 5 0° ± 6° and ext ends only to 
« O'.'OOS from the nucleus (tMarkoff et al.ll2008l '). We plotted 
its direction on Fig[71 Although it points to some blueshifts 
observed beyond ~ l'.'5 to the left of the nucleus in the Fig- 
ure, the gas closer to the nucleus (thus at the scale of the 
jet extent) is not observed in blueshift, suggesting that the 
blueshifts are not related to the radio jet. 

Within the inner arcsecond around the nucleus, the in- 
crease in the gas velocity dispersion may be related to the 
compact nuclear jet at PA = 50° in interaction with the sur- 
rounding gas. Coincidently, this is approximately the PA 
along which we observe an increase in the [01 ] /Hq line ratio 
(Fig.fTOj). a well known signature of shocks l|Baldwin et al.l 
ll98ll : IOsterbrocklll989f) . 



6.3 Estimating the mass accretion flow 

The residual blueshifts observed in Fig.[7]reach high values of 
~ —100 kms~^ in three regions in the far side of the galaxy. 
If we assume that the gas is located in the galaxy plane, this 
means that the gas is inflowing to the nucleus and we can 
thus estimate the rate of mass inflow. In order to do this, 
we assume that the three regions channel gas towards the 
nucleus. Assuming a similar geometry for the three regions, 
the total gas mass inflow will be three times the flow in 
one region. The flux of matter crossing one region can be 
calculated as: 

where Ne is the electron density, v is the streaming velocity 
of the gas towards the nucleus, rup is the mass of the proton, 
r is the cross section radius of the inflow and / is the filling 
factor. For the purpose of the present calculation, we adopt 
a cross section radius of 17 pc and a typical value for the 
filling factor of / = 0.001. 

The velocity residuals in the three blueshifted regions 
on the far side of the galaxy, which we have interpreted 



as due to streaming motions, range from 



)kms 



to 



?a-120kms"^ (Fig.[7ll. As we have assumed that these 
streaming motions occur in the plane of the galaxy, we 
need to correct these values for the inclination of the galaxy 
i = 58.4°. The resulting average velocity for the stream- 
ing motions is i; = — 117kms~^. For A^'e = 500 cm"'' (ob- 
tained from the [SII] ratio) and w = — 117kms~^, we ob- 
tain the value for the total inflow of ionized gas mass of 
<;/)«4.0 X 10"^ M© yr-^ 

We note that the above calculation is only a rough es- 
timate, as it depends on uncertain parameters, such as the 
geometry of the flow. We point out also the fact that the 
estimated inflow rate is in ionized gas and is thus a lower 
limit to the total inflow rate, as it may be only the ion- 
ized "skin" of a more massive inflow in non-emitting neu- 
tral gas. The presence of molecular gas closer than « 300 pc 
fr om the nucleus s e ems t o be ruled out by the observations 
of ISakamoto et all l|200ll ), who found that molecular gas is 
mainly on a "pseudoring" or spiral arm at about 500 pc from 
the nucleus. 

We now compare the estimated inflow rate of ionized gas 
to the mass accretion rate necessary to produce the lumi- 
nosity of the LINER nucleus of M81, calculated as follows: 

Lboi 

where is the efficiency of conversion of the rest mass 
energy of the accreted material into radiation. For LIN- 
ERs it has been concluded that, in most cases, the ac- 
cretion disk is geometrically thi ck, and optically thin 
(jNemmen et all l2006l : lYuanl l2007l ). This type of accre- 
tion flow is known as RIA F (Radiatively Inefficient Ac- 
cretion Flow l|Naravanll2605l ')'). and has a typical value for 
»7«0.01. The nuclear luminosity can be esti mated from the 
X-ray luminosity of Lx =3.2 x 10*'' ergs"^ (|H3[l99i), us- 
ing the approximation that the bolometric luminosity is 
Lb ~ WLx ■ We use these values to derive an accretion rate 
of rn = 5.65 x 10"* M© yr"\ 

Comparing this accretion rate rh with the mass flow of 
ionized gas along the the three regions derived above 0, we 
conclude that <j!> ~ 7 m. Thus the inflow of ionized gas in the 
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inner tens of parsecs of M81 is already of the order of the 
one necessary to feed its AGN. 

6.4 PCA Tomography 

In order to try to extract more information from the com- 
plex gas kinematics of the nuclear region, we have used 
an alternative approach which has been successfully ap- 
plied to revea l smal l scale features around nearby AGN by 
ISteiner et al.l l|2009l ): the method of Principal Component 
Analysis. The method, which is applied directly to the cali- 
brated data cube, allows the separation of the information, 
originally presented in a system of correlated coordinates, 
into a system of uncorrelated coordinates ordered by prin- 
cipal components of decreasing variance. These new coor- 
dinates called "eigenspectra" (eigenvectors which are func- 
tions of wavelength) reveal spatial correlations and anti- 
correlations associated to features that may be associated 
to emission-line regions. The projection of the data onto 
these new coordinates produce images called "tomograms" , 
which represent slices of the data in the space of the eigen- 
vectors. In order to apply this technique, we first corrected 
the data for differential atmospheric refraction and ap- 
plied a Richardson-Lucy deconvolution algorithm (Richard- 
son 1972; Lucy 1974) to the data cube using 6 iterations. A 
larger number of iterations did not improve the results, as 
also verified in pr evious works. A disc ussion about th i s poin t 
can be found in iMorelU et af] (|2010l ) and iMichardI (| 19961 ). 
The real PSF shows a small variation with wavelenght and a 
Gaussian with FWHM = 0.6 " was found to be a good repre- 
sentation of the average PSF over the different wavelenghts. 
We then applied the PCA to the inner 6x14 arcsec^ and, 
in order to enhance the features in the vicinity of the active 
nucleus, also to the inner 4x4 arcsec'^ of the galaxy. We re- 
stricted the spectral range to AA 6200-6700 A, as there is no 
relevant information on the gas emission for A < 6200A 
and the sulphur lines at AA 6717,6731 are affected by falling 
in the gaps of the CCD, as previously explained. 

The resultant eigenspectra and respective tomograms 
did not change significantly when the analysis was per- 
formed on either the inner 6 x 14 arcsec'^ or 4 x 4 arcsec^ of 
the galaxy, the second eigenspectra and tomogram being an 
exception, as they presented small but significant changes. 
Because of the similarity we choose to primarily focus on 
the results corresponding to the larger spatial region for all 
eigenspectra and tomograms except for the second, for which 
both are shown. 

W e have performed a "scree test" (see ISteiner et al.l 
l|2009l) ') to select the significant components - the ones not 
dominated by noise, which we concluded are the first 8 (PC 1 
to PC 8) . The first six eigen-spectra and corresponding to- 
mograms are shown in Figs. 1111 1121 and 1131 PC 7 and PC 8, 
although significant, do not show any particular new struc- 
ture. 

Most of the variance is in the first eigenspectrum PC 1 
(principal component 1, which contributes with 97.50% of 
the variance) , which is similar to the spectrum of the nucleus 
plus central stellar bulge (bottom left panel of Fig[TT} . The 
second tomogram, PC 2 (with 2.13% of the variance), shows 
an anticorrelation between the nucleus and the stellar bulge. 
In the eigenspectrum 2 (bottom central panel of Fig lllf) one 
can clearly see a pure nuclear emission spectrum (positive 



values) anticorrelated with the stellar emission (negative val- 
ues) - also evidenced by the 6500 A stellar absorption fea- 
ture, which appears in absorption in eingenspectrum 1 but 
appears inverted in the eigenspectrum 2. A somewhat con- 
spicuous red Ha wing may be reminiscent of t he double- 
peaked broad Ha reported bv lBower et al.l lll996l). This fea - 
ture is certainly weaker than reported by Bower et al.l (|l996l ) 
and the stellar feature at 6500 A should not be confused with 
emission from the blue Ha wing. When the analysis is per- 
formed on the inner 4x4 arcsec^ (right panels of Fig. lll|l 
a weak feature not seen in the centroid velocity maps ap- 
pears: double peaks in the emission lines of Ha-|-[NII], with 
the blue and red peaks slightly blueshifted («— 80kms~^) 
and slightly redshifted (~80kms~^), respectively, relative 
to the systemic velocity. The fact that the two peaks are 
correlated means that they are spatially coincident, that is, 
they come from the same region and are not spatially re- 
solved. This can be interpreted as due to a compact bipolar 
outfiow or rotation, not resolved in the data. 



6.4-.1 Gas rotation versus outflow: feeding & feedback 

The third, fourth and fifth eigenspectra, PC 3 (with 0.25% of 
the variance), PC 4 (0.07%) and PC 5 (0.01%), show peaks 
associated to the Ha-|-[NII] emission lines spatially anti- 
correlated with symmetric lobes (Fig. ll2p . The pattern in 
the tomogram of PC 5 is, in fact, very similar to that of 
PC 3, and we thus consider that PC 3-f PC 5 correspond to 
one structure and PC 4, whose geometry in the tomogram is 
perpendicular to that of PC 3-l-PC 5, to another structure. 
These patterns can be explained either by a rotating gaseous 
disk or, else, by outflow of ionized gas, either a jet or a milder 
accretion disk outflow. A gaseous disk with such a behavio r 
was already proposed for NGC 4637 bv lSteiner et al. I l|2009t ). 
In the present case, however, we see two such structures 
with axes approximately ortogonal to each other. This can 
be explained by a simple model of a rotating disk and an 
outflow driven by this disk or by an unresolved accretion 
disk inside. This explains not only the existence of the two 
patterns, seen in PC 3-f PC 5 and PC 4, but also the reason 
why they are orthogonal. 

The immediate question that arises is which configu- 
ration, PC 3-l-PC 5 or PC 4, is associated to the rotating 
gaseous disk and which one is produced by the outfiow. We 
propose that PC 3 and PC 5 are associated to the gaseous 
disk and PC 4 to the outfiow. This interpretation is cor- 
roborated by the overall gas rotation seen in the galaxy 
as shown by the [N II] velocity field displayed in Fig.[71 
The rotation axis seen in this figure (close to the nucleus) 
approximately coincides with the one seen in Fig. [12] for 
PC 3-l-PC 5 as does the sense of redshift-blueshift rotation. 
The PC 3 tomogram shows that the disk is within ~ 50 pc 
from the nucleus and comparison of this tomogram with the 
gas kinematics of Fig. [7] suggests that the disk is being fed 
by the gas circulating in the inner part of the galaxy. A 
compact disk in rotation around the nucleus has also been 
found in the H2 kinematics of the inner few tens of parsecs 
of NGC 4151 lIStorchi-Bergmann et al.l l2010t) and Mrk 1066 
(jRiffel fc Storchi-Bergmannll2011^ . 

Why do we see two eingenvectors, PC 3 and PC 5, with 
similar spectral structures? They are not really similar. In 
PC 3 the [O I] emission is seen in correlation with a bluish 
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Figure 11. Left: first eigenspectra (bottom) ad its respective tomograms (top); center and riglit: second eigenspectra and tomogram 
resultant from tlie application of PCA to the larger (central panels) and smaller (right panels) spatial regions. The square in the central 
panel shows the smaller central region where the PCA was also applied. 



continuum while in PC 5 this emission is correlated with a 
reddish continuum. They very probably are associated to 
two regions affected by different dust reddening or, else, by 
distinct scattering properties, or both. Thompson scattering 
by free electrons is wavelength independent but dust scat- 
tering or molecular Rayleigh scattering enhances blue char- 
acteristics. The partially ionized region responsible for the 
[OI] emission seen in tomogram 3 could be exposed to the 
central blue continuum source and scatter it to our line of 
sight, enhancing the blue aspect of the continuum emission. 

We interpret PC 4 as due to an outflow, which most 
probably originates in the interaction of the radio jet dis- 
cussed in Sec. l6.2l with the circumnuclear gas. We plotted the 
direction of the r adio jet on tomogram 4 (arrow in Fie 112|) . 
Considering that iMarkoff et al.l (|2008l ) argues that this jet 
precesses within 20° from this orientation, the jet direction 
approximately coincides with that of the axis of the pre- 
sumed disk in tomogram 3. This is what one would expect 
if a conical outflow or jet is induced by an inner accretion 
disk; assuming that the axis of this inner accretion disk is 
aligned with the structure seen in tomogram 4, we have a 
self-consistent picture. As the highest energy photons are 
probably escaping perpendicular to the disk, we further pre- 
dict that high ionization species, not present in the spectral 
range of our observations, such as [OIII], [Nelll] or Hell 
may be present in the blobs associated to PC 4. 



Tomograms 3 and 4 also show that the [O I] emission 
is spatially correlated with the gas emission to the E-N of 
the nucleus. This is revealed in the eigenspectrum 3 by the 
fact that [O I] is correlated with the blueshifted emission 
(to the east), while in the eigenspectrum 4 it is correlated 
with the redshifted emission (to the north east), in agree- 
ment with the previous result that the [O I] emission is en- 
hanced towards the radio jet orientation of PA ~ 50° . In ad- 
dition, the [O I] profile has indications of double structure, 
being slightly redshifted. Although this question needs fur- 
ther studies we wonder whether this asymmetric [O I] en- 
hancement is not related to the one-sided radio jet, perhaps 
due to shock heating of the partially ionized gas. 

Finally, we show in Fig. [13] the tomogram (left panel) 
and eingenspectrum (right panel) of PC 6. The eigenspec- 
trum shows that the unresolved broad component of Ha 
(see Fig.[8)l, which is confined to the red circle at the nu- 
cleus in the tomogram, is spatially anti-correlated with the 
narrow-emission lines, which originate in the blue patches 
surrounding the nucleus. 

7 CONCLUSIONS 

We have measured the kinematics of the stars and ionised 
gas in the inner 120 x 250 pc^ of M81 using integral field 
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Figure 13. Left: Tomogram of PC 6; right: eigenvector 6, showing the broad Ha component. 



spectroscopy obtained with the GMOS instrument at the 
Gemini North telescope with a spatial resolution of 10 pc at 
the galaxy. The main results of this paper are: 

• The stellar velocity field (centroid velocity) shows circu- 
lar rotation with a major axis orientation of ~ 153° , which is 



similar to that of the large scale stellar rotation (as reported 
in previous works). But the two-dimensional coverage al- 
lowed the detection of deviations from pure rotation which 
can be attrib uted to stellar motions possibly associated to 
a nuclear bar (jElmegreen et al.ll 19951 ): 

• The stellar velocity dispersion of the bulge is 
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162 ± 15 km s ^ , in good agreement with previous measure- 
ments, resulting in an estimate of the black hole mass of 
Mbh = 5.5l^;^ xlO'^Mo; 

• The gaseous velocity field is completely distinct from 
the stellar one, and is dominated by blueshifts in the far 
side of the galaxy and a few redshifts in the near side. The 
subtraction of the stellar velocity field from the gaseous one 
confirms the presence of excess blueshifts of w — lOOkms"^ 
on the far side and some redshifts in the near side. If the gas 
is in the galaxy plane, these excesses can be interpreted as 
streaming motions towards the center of the galax;y; 

• Prom the measured gas velocities and assumed ge- 
ometry, we estimate an ionized gas mass flow rate of 

4.0 X 10^"^ M© year^^, a value about seven times the 
necessary accretion rate to power the AGN. However, this 
is only the mass flow rate in ionized gas. Most probably, the 
total inflow is much larger but is dominated by non-emitting 
gas; 

• We have applied the technique of principal component 
analysis (PCA) to our datacube and find that the PGA re- 
veals small scale features not seen in the measured velocity 
field. The main features are: (1) an unresolved bipolar out- 
flow or rotation at the nucleus; (2) a nuclear rotating disk 
within « 50 pc from the nucleus; (3) a compact bipolar out- 
flow approximately perpendicular to the disk; 

• The gas velocity field suggests that the nuclear disk is 
being fed by gas circulating in the galaxy plane; 

• An increase in the [01] /Ha line ratio occurs at approx- 
imately the same orientation of a one-sided nuclear radio jet 
(PAw50°), which is also the approximate orientation of the 
bipolar outflow revealed by the PGA. The increase in the 
[01] /Ha line ratio is thus probably due to the interaction 
of the jet with the circumnuclear gas, as this ratio is a well 
known signature of shocks. 

With our observations we are thus resolving both the 
feeding of the M 81 AGN - via the compact disk and sur- 
rounding gas which seems to be replenishing the disk, and 
its feedback - via the compact outflow approximately per- 
pendicular to the disk, probably driven by this disk or by 
an unresolved accretion disk inside. 
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